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ABSTRACT 

The  purpose  of  this  project  is  to  localize  buried 
objects  in  water-saturated  sand  by  variable  incidence 
acoustic  pulse  reflections.   In  particular,  thin  aluminum 
plates  are  used  to  model  the  buried  objects.   The  sound 

speed  and  the  depth  of  the  overlying  sand  layer  are  predicted 

2   2 

by  using  the  T  -X   and  the  ray  parameter  methods.   The 

experimentally  computed  results  have  an  average  3.5%  depth 

error.   The  calculated  speed  of  sound  in  the  water-saturated 

2   2 
sand  layer  differed  by  8.1%  between  the  T  -X   and  the  ray 

parameter  methods . 
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I.   INTRODUCTION 

A.   OBJECTIVES 

Identification  of  buried  objects  in  the  ocean  bottom  is 
an  important  problem  in  many  fields  such  as  oil  exploration, 
geophysical  research,  and  ocean  engineering.   A  simpler 
problem  is  the  identification  of  layers  of  infinite  extent 
and  unknown  properties  in  the  ocean  sediments.   The  objec- 
tives of  this  thesis  are: 

1.  to  determine  experimentally  the  depth  of  burial 
of  a  thin  plate  in  a  simulated  sedimentary  bottom; 

2.  to  determine  the  sound  speed  and  the  thickness 
of  the  overlying  sand  layer;  and, 

3.  to  determine  the  sound  speed  and  thickness  of 
the  thin  plate. 

It  was  the  intent  of  this  thesis  to  determine  all  of  these 
properties  by  remote  sensing  with  acoustic  pulses  in  order 
to  simulate  an  ocean  measurement  with  a  ship  and  near- 
surface  sensors. 

The  first  objective  could  have  been  accomplished  by 
normal-incidence  time-of-f light  measurements;  however,  the 
second  objective  requires  more  information  than  that  so  a 
variable  angle  of  incidence  technique  was  used.   This  second 
objective  was  actually  the  principal  goal  since  the  third 
objective  could  not  be  achieved  by  using  the  technique 
described  herein.   The  technique  failed  because  the  pulse 
could  not  be  made  short  enough  to  distinguish  the  reflection 


from  the  plate's  lower  surface.   This  is  primarily  the 
result  of  using  very  thin  plates  with  respect  to  the  acoustic 
wavelength.   Perhaps  the  amplitude  of  the  reflection  from 
the  top  of  the  plate  could  be  used  to  assist  in  identifying 
the  plate  material  but  the  overlying  sand  layer  makes  this 
very  difficult. 

B .  APPROACH 

The  approach  used  is  as  follows: 

1.  Short  pulses  consisting  of  several  cycles  of  a 
150kHz  sinusoidal  wave  were  used  as  the  source 
pulse. 

2.  A  directional  projector  was  used  to  minimize  inter- 
ference from  surface  and  side-wall  reflections. 

3.  The  receiving  hydrophone  was  mounted  on  a  carriage 
so  that  the  horizontal  distance  from  source  to 
receiver  could  be  easily  varied . 

4.  Each  of  the  thin  plates  were  buried  in  the  sand 
bottom  of  a  small  water  tank. 

5.  The  plates  were  aluminum  ranging  in  thickness  from 
0.1  to  1.8  cm  and  were  buried  at  depths  from  1.8 
to  5.0  cm. 

6.  Two  analytical  methods  were  used  for  reducing  the 
data.   One  method  was  based  on  the  relationship 
between  the  square  of  the  travel  time  (T)  of  the 
reflections  and  the  square  of  the  source-to-receiver 
horizontal  range  (X) .   The  other  method  was  based 

on  maintaining  a  constant  ray  parameter  (P)  defined 
as  the  sine  of  the  incident  angle  divided  by  the 
local  sound  speed. 

C.  BACKGROUND 

Reference  1  gives  an  excellent  historical  development 
of  various  seismological  methods  for  determining  properties 
of  sediment  layers.   The  following  is  a  brief  summary  of 
that  account. 
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Karcher  and  Haseman  made  the  first  successful  seismic 
reflection  measurements  in  1921.   Techniques  for  acoustic 
determination  of  sediment  properties  have  advanced  greatly 
since  then  and  generally  fall  into  one  of  three  categories: 
normal  incidence,  wide-angle  reflection,  and  refraction 
methods.   Most  of  the  development  has  been  done  by  the 
petroleum  exploration  industry;  however,  the  geophysical 
community  has  used  these  methods  to  study  the  earth's  crust 
as  early  as  the  late  1930's.   Dix,  Krey  and  Durbaum 
refined  the  sound  speed  measurements  for  sediment  layers 
below  the  first  layer  by  properly  accounting  for  the  shape 
of  the  arrival  time  as  a  function  of  range  curves .   By 

19  55,  Dix  had  demonstrated  the  equivalence  of  T-X  methc 

2   2 

and  T  -X  method  and,  in  196  8,  Le  Pichon  refined  the 

2   2 
T  -X  method  until  it  became  a  standard  analysis  tool. 

A  major  enhancement  was  introduced  by  Bryan  (see  Ref .  2)  in 

2   2 
19  74  with  the  ray  parameter  method.   Both  the  T  -X  method 

and  the  ray  parameter  method  were  used  in  the  data  analysis 

for  this  thesis. 
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II.   THEORY 

2   2 

A.   T-X  AND  T  -X   METHODS 

A  simple  model  of  a  layered  sediment  is  shown  in  Fig. 

2.1.   The  layers  have  thicknesses  H, ,  H      and  H-,  and  sound 

2   2 

speeds  V1 ,  V-  and  V3 .   The  T  -X   method  [Ref.  1]  can  be 

applied  to  this  situation  as  described  in  this  section. 


Figure  2.1.   Geometry  of  Layers 

For  convenience,  let  the  sound  speed  in  the  uppermost 
layer  (the  water  layer)  be  less  than  that  in  the  first 
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sediment  layer.   The  length  of  the  path  that  reflects  from 
the  top  of  the  sediment  layer  is  given  by, 


Rll   =   x2  +  4<Hi  "s>2  <!) 


To  convert  this  to  travel  time  squared,  divide  through 
by  V2: 

T2(X)   =   T2(0)  +  (ix2  (2) 

Vl 

where  T(X)  is  the  reflection  travel  time  at  horizontal  range 

X  and  T(0)  is  the  intercept  of  the  travel  time  curve  at 

2 
X  =  0.   Since  T  (0)  is  a  constant,  the  reflection  curve  is 

a  hyperbola.   The  corresponding  curve  for  the  direct  path 

from  the  source  to  the  receiver  is  a  straight  line  and  this 

hyperbola  is  asymptotic  to  that  direct-path  line. 

2   2 
When  plotted  in  the  T  -X   plane,  however,  these  lines 

are  both  straight.   In  this  way,  the  slope  of  the  reflection 

curves  can  be  easily  determined.   As  can  be  seen  from  the 

2   2 
equation  for  the  reflection  path,  the  slope  in  the  T  -X 

plane  is  — ~-  which,  in  this  thesis  work,  is  calculated  by 

Vl  2   2 

linear  regression  of  the  T  -X   data.   This  rather  simple 

but  accurate  method  is  due  to  Green  [Ref .  1]  and  is  asso- 
ciated with  the  wide-angle  seismic  reflection  technique. 

The  paths  that  pass  through  one  or  more  sediment  layers 
before  reflection  (see,  for  example,  R, „  in  Fig.  2.1)  do  not 
lead  exactly  to  hyperbolas  but  the  departure  is  not  great 
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A. 


B. 


^V-C*) 


Figure  2.2.   T-X  and  T2-X2  Curves 
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and  only  a  small  amount  of  error  is  introduced  by  fitting 

2   2 
this  data  with  straight  lines  in  the  T  -X   plane. 

If  the  layers  were  extremely  thick  or  had  great  differ- 
ences in  sound  speed,  then  the  error  might  become  apprecia- 
ble.  The  reciprocal  of  the  slope  of  the  R,  2  curve  is  equal 
to  the  average  speed  (squared)  through  both  the  water  and 
the  first  sediment  layer: 


V   =   (H1+H2)/((H1/V1)  +  (H2/V2)) 


=   (V1T1+V2T2)/(T1+T2)  (3) 

Representative  T-X  curves  for  this  case  are  shown  in 
Fig.  2. 2. a.   Curve  D  represents  the  direct  path  and  R, , 
represents  the  reflection  curve  for  reflection  from  the  top 
of  the  sediment.   For  R.  ~ ,  the  exact  curve  (solid  line)  is 

given  as  well  as  the  hyperbola  approximation  (dotted  curve) 

2   2 
to  that  curve.   The  equivalent  curves  in  T  -X   space  are 

shown  in  Fig.  2.2.b.   Here  the  inverse  slopes  give  the  sound 

speed  squared.   The  intercept  T(0)  gives  the  two-way 

vertical  travel  time  through  the  layers. 

2   2 
The  principal  limitation  of  the  T  -X  method  is  shown 

in  these  figures.   For  any  reflection  other  than  the  reflec- 
tion from  the  top  of  the  sediment,  the  curve  is  not 

2   2 
exactly  a  straight  line  in  T  -X   coordinates  but  it  falls 

slightly  below  the  equivalent  straight  line  as  range 

increases . 
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B.   THE  RAY  PARAMETER  METHOD 

In  order  to  measure  the  properties  of  a  layer  without 
knowing  the  structure  of  the  overlying  layers,  the  ray 
parameter  method  can  be  used  [Ref .  2] .   Any  ray  path  is 
specified  by  a  fixed  value  of  the  ray  parameter  P  defined 
by, 

P   =   sin  9  (z)/c(z)  (4) 

where  z  is  depth,  c(z)  is  the  sound  speed  profile  and  6(z) 
is  the  angle  that  the  ray  makes  with  the  vertical  at  depth 
z.  Figure  2.3  shows  the  relationship  between  the  incident 
angle  and  the  ray  parameter. 

According  to  the  ray  parameter  method,  the  increment  of 
travel  time  along  the  ray  path  gives  the  corresponding 
increments  in  depth  and  range: 


2  2  1/2 
dz   =   c  dTcos  6   =   (1  -P  c  )    c  dT  (5) 


—  1  /  9 

dX   =   c  dTcos  9   =   P  c2dt   =   Pc(l  -Pc)     dz      (6) 


After  integrating  over  depth  from  the  surface  to  depth 
z, ,  these  expressions  can  be  rewritten  as  follows, 


T(P)   =    /   c"1(l  -  P2c2)"1/2  dz  (7) 

0 
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Figure  2.3.   Geometry  of  Ray  Parameter 


r1  2    2    ~1/2 

X(P)   =    /    Pc(l  -Pc)     dz 

0 


(8 


These  equations  give  the  travel  time  T  and  the  hori- 
zontal distance  X  along  the  ray  path  given  by  P  in  terms  of 
integrals  over  depth.   The  equations  must  be  solved  for  the 
sound  speed  variation  c(z).   This  can  be  done  for  some 
special  cases  of  layer  structure. 

Consider  a  homogeneous  layer  (sound  speed  c, )  extending 
from  z,  to  z,  +  Az  (Fig.  2.3).   For  the  ray  of  a  particular 
ray  parameter  value,  there  are  two  (in  this  case)  rays  that 
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return  from  the  bottom  with  that  same  parameter  value 
although  they  arrive  at  the  receiver  at  different  ranges. 
Note  that  these  rays  are  parallel  through  the  water  column. 
The  travel  time  and  horizontal  distance  for  these  two  rays 
are,  according  to  Eqs .  (7)  and  (8), 


Zl+Az   _  -1/2 

AT(P)   =   2    /       c"1(l-P2c2)     dz  (9) 

Zl 


z1+Az  _1/2 

AX(P)   =   2    /       Pc(l  -P2c2)     dz  (10) 

Zl 


where  the  factor  2  accounts  for  the  round  trip  of  the  ray. 
If  the  sound  speed  is  constant  at  c,  in  the  layer,  these 
equations  can  be  integrated  to  give, 


-1      2  2  "1/2 
AT   =   2c,  (1  -Pc)     Az  (11) 


2  2  "1/2 
AX   =   2Pc1(l  -P  cp     Az  (12) 


The  expressions  for  sound  speed  and  layer  thickness  are 
then, 


cl  =   ^AT*172  (13) 


1  2  2  1/2 

Az   =   |  c,AT(l  -P  cp  (14) 
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Figure  2.4.   T-X  Curve  and  the  Ray  Parameter 

The  important  feature  of  the  ray  parameter  method  is  that 
the  structure  above  the  layer  has  dropped  out  because  of 
the  choice  of  ray  path. 

In  order  to  apply  this  technique  in  practice,  an  initial 
value  of  ray  parameter  was  chosen  based  on  a  geometry  that 
gives  good  clean  reflection  signals.   After  the  T-X  curves 
are  drawn,  this  value  of  ray  parameter,  which  gives  the 
slope  of  the  T-X  curves  at  the  points  of  interest,  was  used 
to  identify  the  two  parallel  return  rays.   The  principal 
difficulty  of  this  method  is  that  the  T-X  curves  are  not 
straight  lines  and,  therefore,  curves  must  be  fit  to  the 
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experimental  data.   This  introduces  some  uncertainty  into 
the  process  but,  fortunately,  the  values  for  layer  thickness 
is  not  very  sensitive  to  the  type  of  fit  that  is  done. 
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III.   EXPERIMENTAL  PROCEDURES 

A.   ARRANGEMENT  OF  THE  EXPERIMENT 

The  measurements  were  made  in  a  tank  which  was  300  cm 
long,  115  cm  wide  and  90  cm  deep.   The  tank  was  filled  with 
carefully  de-aerated  (fresh)  water  and  had  a  55  cm  deep 
sand  bottom  (also  de-aerated) . 

At  first,  an  ITC  501  projector  was  used  but  the  beam 
pattern  was  fairly  broad  leading  to  excessive  reflections 
from  the  surface  and  the  side-walls  and,  more  importantly, 
this  projector  was  large  enough  so  that  it  would  ring  for  a 
significant  length  of  time  compared  to  the  desired  pulse 
length.   As  a  result,  the  output  pulse  was  too  long  regard- 
less of  the  length  of  the  pulse  used  as  input. 

These  problems  led  to  the  selection  of  a  USRF  F-41 
projector,  which  is  a  directional  laboratory  standard  trans- 
ducer used  in  acoustic  reciprocity  calibrations  over  the 
frequency  range  from  15  to  150  kHz.   This  projector  consists 
of  lead  zirconate-titanate  elements  cemented  to  a  high- 
density  kennametal  disks.   The  resulting  array  of  elements 
is  about  3.8  cm  wide  and  5  cm  high. 

The  receiving  hydrophone  used  was  an  omnidirectional 
LC-10  hydrophone  mounted  on  a  carriage  that  was  free  to  roll 
along  tracks  installed  on  the  long  sides  of  the  tank.   In 
this  way  the  horizontal  distance  (X)  could  be  easily  varied. 
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A  block  diagram  of  the  electronics  used  is  shown  in 
Fig.  3ol.   The  digital  oscilloscope  was  used  to  measure  the 
time  differences  between  various  pulses.   In  order  to  reject 
unwanted  noise,  a  filter  was  used  to  pass  only  100  kHz  to 
200  kHz. 

The  choice  of  pulse  length  was  governed  by  the  criteria 
that  it  be  (1)  long  enough  to  reach  steady-state  conditions 
in  the  projector  and  (2)  short  enough  to  resolve  the  indi- 
vidual reflections  from  the  layers.   The  best  compromise 
turned  out  to  be  four  cycles  of  a  150  kHz  sine  wave,  but, 
unfortunately,  this  was  not  short  enough  to  completely 
define  the  thin  plate.   At  150  kHz,  one  cycle  is  6.7  yS 
so  the  four-cycle  pulse  was  about  27  yS  long.   A  typical 
received  signal  using  the  four-cycle  pulse  is  shown  in 
Fig.  3.2  while  Fig.  3.3  shows  the  same  situation  with  an 
eight-cycle  pulse.   The  four-cycle  pulse  gives  much  better 
definition  of  individual  arrivals. 

Before  the  actual  measurements  were  made,  several  tests 
were  performed  in  order  to  insure  that  side-wall  and  surface 
reflections  would  not  interfere  with  the  measurements.   The 
depths  of  the  projector  and  hydrophone  were  determined  by 
trial  and  error  not  only  to  minimize  interfering  reflections 
but  also  to  optimize  the  separation  between  the  desired 
arrivals.   In  addition,  this  set  a  practical  limit  on  the 
horizontal  range. 

An  attempt  was  made  to  move  the  axis  of  the  source-to- 
receiver  line  toward  the  tank  wall  so  that  two  parallel 
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Figure  3.3.   8  Cycle  Pulse  Signal 

measurements  could  be  made  without  disturbing  the  sand  to 
bury  each  plate.   This  introduced  unacceptable  side-wall 
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Figure  3.5.   Projector  Is  Not  on  the  Center  Line 
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reflections  so,  instead,  three  plates  were  buried  in  line 
along  the  long  axis  of  the  tank  with  the  projector  and 
receiver  positioned  as  shown  in  Fig.  3.4. 

Here,  the  plates  were  buried  at  different  depths  and  the 
projector  and  hydrophone  were  positioned  as  shown.   (Notice 
that  in  the  third  set-up  on  the  right  the  source  and  receiver 
are  reversed  to  avoid  reflections  from  the  end-wall.) 

Figure  3.5  shows  the  side-wall  reflections  that  are 
introduced  when  the  source-to-receiver  axis  is  displaced 
only  15  cm  toward  one  side  of  the  tank.   This  result  should 
be  comparable  to  Fig.  3.2  but  only  a  position  along  the 
centerline  of  the  tank  turned  out  to  be  acceptable. 

B.   EXPERIMENTAL  PROCEDURE 

Once  the  experiment  was  set  up,  the  first  step  was  to 
identify  the  individual  arrivals  from  a  particular  buried 
plate.   By  considering  the  projector  depth  S,  the  hydrophone 
depth  h,  the  water  depth  H,  the  major  incident  angle,  and 
the  slant  distance  R,  the  time  of  arrival  of  the  various 
reflections  was  determined  by  simple  image  theory  [Ref.  3, 
pp.  427-429].   A  representative  diagram  of  the  various 
arrivals  is  shown  in  Fig.  3.6  for  constant  sound  speed  in 
the  water  layer. 

Since  a  directional  projector  was  used  and  was  aimed 
downward  toward  the  sediment,  the  surface-reflected  path 
7  could  be  neglected.   For  most  of  the  measurements,  paths  1, 
3  and  4  were  used. 
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Figure  3.6.   Waves  in  the  Isospeed  Shallow 
Water  Channel 


Basically,  three  methods  were  used  to  initially  identify 
the  arrivals  on  the  digital  oscilloscope.   First,  the 
arrival  times  predicted  by  the  method  of  images  were  com- 
pared to  the  measured  arrival  times.   The  amplitude  of 
the  arrivals  also  helps  in  the  identification  process  since 
the  pressure  amplitude  is  reduced  by  the  reflection  coeffi- 
cient at  each  of  the  interfaces  of  the  layers.   Finally,  a 
path  that  reflects  from  the  surface  can  be  readily  identi- 
fied by  disturbing  the  water's  surface  with  a  ruler.   Any 
surface-reflected  arrivals  change  greatly  on  the  oscillo- 
scope display  and  can  easily  be  separated  from  those 
arrivals  that  only  reflect  from  the  bottom. 
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A  typical  display  on  the  oscilloscope  is  shown  in  Fig. 
3.7.   Here,  the  projector  is  inclined  so  that  the  beam  has 
an  incident  angle  of  55  degrees  and  the  water  depth  is  35  cm 
The  arrivals  shown  are: 

1.  direct  wave, 

2.  surface-reflected  wave, 

3.  bottom-reflected  wave, 

4.  wave  reflected  from  the  top  of  the  plate, 

5.  wave  reflected  from  the  bottom  of  the  plate  after 
passing  through  the  plate, 

6.  bottom-surface-reflected  wave. 


1      2 


3  4    5     6 


Figure  3.7.   Different  Kinds  of  Waves  on  Oscilloscope 
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Once  the  arrivals  are  identified  for  one  source-to-receiver 
range,  the  signals  can  be  followed  by  eye  if  the  receiver 
carriage  is  moved  smoothly. 

Next,  the  range  X  is  varied  by  moving  the  carriage  and 
the  arrival  times  T  of  the  individual  reflections  are 

recorded.   These  measurements  form  the  basis  of  the  T-X  and 

2   2 
T  -X  methods.   When  these  measurements  are  plotted  as 

arrival  time  T  as  a  function  of  range  X,  the  direct  arrivals 

appear  as  straight  lines  and  the  reflected  arrivals  (either 

from  the  top  of  the  sand  layer  or  from  any  of  the  layers 

below)  appear  as  parabolas.   (They  are  actually  closer  to 

hyperbolas  but,  over  the  limited  range  available  in  this 

experiment,  parabolas  were  used  as  an  approximation.) 

2 
The  T  and  X  measurements  were  also  plotted  as  T   as  a 

2 
function  of  X  .   On  these  coordinates,  any  of  the  lines  are 

straight  (or  nearly  so).   From  the  slope  of  these  lines, 

the  sound  speed  either  in  the  water  or  in  the  layers  could 

be  calculated.   A  linear  regression  method  was  used  to 

compute  the  slope  for  the  individual  arrivals. 

Next,  the  ray  parameter  P  was  calculated.   This  value 

is  fixed  for  a  given  measurement.   The  incident  angle  was 

determined  from  the  source-receiver  geometry  and  the  T-X 

data  were  fit  with  parabolas.   Then  P  could  be  calculated 

2 

directly  from  the  parabola  equation:   T  =  aX   +  bX  +  c , 

dT/dX  =  P  =  2aX  +  b. 

Once  the  parabolas  have  been  determined,  the  layer  depth 
can  be  found  from  the  T-X  measurements.   For  example,  the 


29 


arrival  from  the  wave  that  passed  through  the  upper  sand 

layer  would  be  described  by  a  parabola  of  the  form 

2 
Tl  =  alXl  +  blXl  +  cl  and  the  reflection  from  the  top  of  the 

2 

sand  layer  would  be  given  by  T2  =  a2X   +  b2X   +  c    .      The 

ray  parameter  P  must  be  the  same  for  each  arrival;  there- 
fore, for  the  wave  passing  through  the  sand  layer, 
P  =  2a,  X..  +  b,  or  X..  =  (P-b,)/2a, ,  and,  for  the  wave  reflected 
from  the  top  of  the  sand  layer,  P  =  2a2X2  +  b2  or  X2  =  (P-b9)/2a2. 
Since  AX  =  X,-X9, 


P-b,  2      P-b,  P-b9  9      P-b9 

AT    =     [Cai(^7»     +  bi(T57>   +  ci}  "  {a2(-2?7)    +b2(TiT»+c2}] 


These  two  values  (T  and  X)  are  necessary  to  calculate  the 

2   2 

depth  and  the  sound  speed  (to  compare  with  the  T  -X  deter- 
mination of  sound  speed)  of  the  layer.  The  values  of  layer 
depth  Az  and  sound  speed  c,  are  given  by, 

1  2  2  1/2 

Az   =   |  c,AT(l  -P  cj) 


2   2 
The  first  value  of  c.  is  from  the  T  -X  method  and  AT  is 

from  the  parabolic  fit  of  the  T-X  data.   The  new  value  of 

c,  can  be  compared  to  the  old  value  as  a  check  on  the  two 

methods . 

C.   SUMMARY  OF  EXPERIMENTAL  PROCEDURE 

1.  Bury  the  thin  homogeneous  plates  at  various  depths; 

2.  Identify  the  arrivals  on  the  oscilloscope; 
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3.  Measure  the  horizontal  distance  X  and  the  arrival 
time  T  of  the  various  arrivals; 

2   2 

4.  Make  the  T-X  and  T  -X  plots; 

5.  Find  the  sound  speed  c,  in  the  layer  by  using  the 
T2-X2  method; 

6.  Determine  the  incident  angle  from  the  source- 
receiver  geometry; 

7.  Calculate  the  ray  parameter  P  (fixed  for  each 
experiment) ; 

8.  Fit  the  T-X  data  with  parabolas; 

9.  Calculate  AX  and  AT  by  the  ray  parameter  method; 

10.  Calculate  the  layer  depth  Az; 

11.  Calculate  c,  and  compare  to  the  value  determined 
previously. 
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IV.   ANALYSIS  AND  RESULTS  OF  EXPERIMENTS 

A.   MEASUREMENT  AND  ANALYSIS 

A  measurement  of  the  sound  speed  in  the  water  was 

required  for  each  of  the  reflection  analyses  so  this  was 

2   2 

measured  using  the  T  -X  method.   The  slope  of  the  straight 

line  fit  to  the  data  points  gave  a  sound  speed  of  1447  m/s 
and  the  correlation  coefficient  associated  with  the  linear 
regression  was  almost  100%.   For  this  measurement,  the  pro- 
jector beam  was  aimed  horizontally  and  the  hydrophone  was 
positioned  at  the  same  depth  as  the  projector. 

In  each  of  the  subsequent  experiments,  the  projector  and 
hydrophone  depths  were  chosen  to  optimize  the  separation 
and  clarity  of  the  reflection  signals.   In  every  case,  the 
source  signal  was  four  cycles  of  a  150  kHz  sine  wave. 

1.   Measurement  1 

For  this  experiment  an  aluminum  plate  (thickness  of 
1.7  cm)  was  buried  3.5  ±  0.05  cm  below  the  sand's  surface 
and  the  projector  depth  was  8.5  cm  while  the  hydrophone 
depth  was  6.5  cm. 

TABLE  I 
Speed  of  Sound  (Measurement  1) 

WAVE  SPEED  OF  SOUND  (m/s )     CORRELATION  (%)     ERROR  (%) 
SAND             1661                  98.3 

BOTTOM  1415  99.4  2.2 

DIRECT  1463  99.8  1.1 
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Table  I  gives  the  speed  of  sound  calculated  from  the  various 
arrivals  along  with  the  correlation  coefficients  of  the 
linear  regressions.   The  speed  of  sound  in  the  water  should 

be  compared  to  the  value  of  14  4  7  m/s  discussed  above.   The 

2   2 
T-X  and  T  -X   data  are  shown  in  Fig.  4.1. 

As  discussed  in  the  previous  chapter,  the  T-X  data 

were  fit  with  a  parabola  for  simplicity.   The  data  for  the 

wave  that  reflected  from  the  top  of  the  aluminum  plate  was 

fit  with 


T   =   0.0926  X2  -  2.89  X  +  459 


and  the  data  for  the  wave  that  reflected  directly  from  the 
top  of  the  sand  was  fit  with 


T   =   0.113  X2  -  3.03  X  +  410 


The  geometry  selected  for  the  ray  parameter  determination 

gave  an  incident  angle  of  29  degrees  which  corresponds  to  a 

—  fi 
ray  parameter  P  equal  to  3.32  x 10    s/m.   This  value  of  slope 

locates  a  point  on  each  of  the  two  parabolas.   The  differ- 
ence in  X  values  gives  a  AX  of  5.47  cm.   The  parabolic  form 
for  the  reflection  from  the  top  of  the  plate  gives  the  time 
difference  AT  of  51.4  yS. 

Next,  the  layer  thickness  Az  was  calculated  from 
Equation  (14).   In  this  case,  Az  was  determined  to  be  3.56  cm, 
This  differs  by  only  0.3%  from  the  measured  value  for  the 
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Figure  4.1.   T-X  and  T  -X   Plots  (Measurement  1) 
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overlying  sand  layer.   (The  burial  depth  of  the  plate  gives 
the  thickness  of  the  top  sand  layer.)   The  value  for  the 
layer's  sound  speed  was  also  calculated  from  the  ray  param- 
eter method  but  its  value  was  not  significantly  different 

2   2 

from  that  calculated  by  the  T  -X  method. 

2 .   Measurement  2 

For  the  second  measurement,  the  hydrophone  depth  was 
6.4  cm,  the  projector  depth  was  8.5  cm  and  the  aluminum  plate 
(1.7  cm  thick)  was  buried  at  a  depth  of  2.9  ±  0.1  cm  in  the 

sand.   The  calculated  sound  speeds  are  given  in  Table  II  and 

2   2 
the  T-X  and  T  -X  data  are  shown  in  Fig.  4.2. 


TABLE  II 
Speed  of  Sound  (Measurement  2) 


WAVE 

SPEED 

OF  SOUND 

(m/s) 

CORRELATION 

(%) 

ERROR  ( 

SAND 

1647 

94.3 

BOTTOM 

1419 

99.6 

1.96 

DIRECT 

1437 

99.1 

0.69 

Figures  4.3  and  4.4  show  representative  displays  from 
the  oscilloscope  at  two  of  the  ranges  at  which  measurements 
were  taken. 

The  direct  wave  (A) ,  the  surface-reflected  wave  (B) , 
the  wave  that  reflects  from  the  top  of  the  sand  (C) ,  the 
wave  that  reflects  from  the  top  of  the  buried  plate  (D) ,  and 
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Figure  4.2.   T-X  and  T  -X   Plots  (Measurement  2) 
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Figure  4.3.   Measurement  2  at  X  =  34.5  cm, 
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the  wave  that  reflects  from  the  top  of  the  sand  and  again 
from  the  water  surface  (E)  are  shown  in  these  figures.   The 
vertical  line  marks  the  beginning  of  the  wave  that  is 
reflected  from  the  top  of  the  plate  and  the  number  gives  the 
arrival  time. 

In  this  case,  the  curve  for  the  arrival  reflected 
from  the  top  of  the  plate  is, 


T   =   0.473  X2  -  32.0  X  +  1020 


and  the  curve  for  the  arrival  reflected  from  the  top  of  the 
sand  is, 


T   =   0.136  X2  -  5.72  X  +  471 


—  f> 
Here,  the  incident  angle  was  28.3  degrees  so  P  was  3.30  x 10 

s/m.   The  ray  parameter  calculations  gave  3.87  cm  for  AX  and 

54  yS  for  dT  so  the  layer  depth  Az  was  3.73  cm.   The  actual 

sand  layer  thickness  was  3.4  cm  so  the  error  was  slightly 

less  than  10%.   The  value  for  sound  speed  in  the  water  by  the 

ray  parameter  method  was  1474  m/s  which  is  about  10%  high. 

3 .   Measurement  3 

In  this  measurement,  an  aluminum  plate  0.1  cm  thick 

was  buried  to  a  depth  of  2.3  ±  0.1  cm.   The  projector  depth 

was  8.5  cm  and  the  hydrophone  depth  was  5.1  cm.   The  result- 

2   2 
ing  speed  of  sound  measurements  (determined  by  the  T  -X 

method)  are  summarized  in  Table  III  and  the  T-X  data  and  the 

38 


2   2 

T  -X  data  are  shown  in  Fig.  4.5.   Figures  4.6  and  4.7 

show  representative  displays  from  the  oscilloscope  at  two 
of  the  ranges  at  which  measurements  were  taken. 

TABLE  III 
Speed  of  Sound  (Measurement  3) 

WAVE  SPEED  OF  SOUND  (m/s)    CORRELATION  (%)    ERROR  (%) 
SAND          1649                    95.2 

BOTTOM        1471  99.9  1.7 

DIRECT        1510  99.8  4.3 

The  parabola  for  the  arrivals  reflecting  from  the 
top  of  the  plate  was  determined  to  be 


T   =   0.290  X2  -  13.3  X  +  580 


and  the  parabola  for  the  arrivals  directly  reflected  from 
the  top  of  the  sand  is 


T   =   0.0484  X2  +  0.528  X  +  357 


For  the  optimum  geoemtry,  the  incident  angle  was  27  degrees 

—  fi 

so  that  the  ray  parameter  P  was  3.00  xio 

The  ray  parameter  analysis  gave  2.53  cm  for  AX  and 
33.38  ys  for  AT;  therefore,  the  sand  layer  depth  was  computed 
to  be  2.42  cm:   less  than  one  percent  from  the  actual  value. 
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Figure    4.5. 


2      2 

T-X    and    T    -X 


Plots  (Measurement  3) 
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Figure  4.6.   Measurement  3  at  X  =  21.7  cm, 


Figure  4.7.   Measurement  3  at  X  =  29.3  cm 
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The  calculated  value  for  the  sound  speed  in  the  sand  layer 

was  1578  m/s  and  this  differs  by  only  3  percent  from  the 

2   2 
value  obtained  by  the  T  -X  method. 

4 .   Measurement  4 

For  this  measurement,  an  aluminum  plate  0.3  cm  thick 

was  buried  to  a  depth  of  1.8  ±  0.1  cm.   The  projector  depth 

was  8.5  cm  and  the  hydrophone  depth  was  3.5  cm.   The  speed 

of  sound  measurements  is  summarized  in  Table  IV  and  the  T-X 

2   2 
data  and  the  T  -X   data  are  shown  in  Fig.  4.8.   Representative 

displays  from  the  oscilloscope  at  two  ranges  are  shown  in 

Figures  4.9  and  4.10. 


TABLE  IV 
Speed  of  Sound  (Measurement  4) 


WAVE 

SPEED  OF  SOUND 

(m/s) 

CORRELATION 

(%) 

ERROR  (%) 

SAND 

1582 

99.7 

BOTTOM 

1458 

99.7 

0.7 

DIRECT 

1461 

99.9 

1.0 

The  parabolic  form  for  the  reflections  from  the  top 
of  the  plate  is 

T   =   0.0259  X2  +  0.361  X  +  429 

and  the  form  for  the  reflections  from  the  top  of  the  sand  is 

T  =  0.0328  X2  +  1.17  X  +  380 
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Figure    4.8.       T-X    and 


2       2 

T-X 


Plots  (Measurement  4) 
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Figure   4.9.      Measurement   4    at   X   =    35.1    cm. 
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Figure    4.10.      Measurement    4    at   X   =    47.7   cm. 
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For  the  optimum  geometry,  the  incident  angle  was  50  degrees 
so  that  the  ray  parameter  P  was  5.3  x 10 

The  ray  parameter  analysis  gave  5.80  cm  for  AX  and 
4  0.4  ys  for  AT;  therefore,  the  sand  layer  depth  was  computed 
to  be  1.74  cm.   The  calculated  value  for  the  sound  speed  in 
the  sand  layer  was  1647  m/s .   In  both  cases,  the  results 

compare  very  well  with  the  measured  value  for  layer  depth 

2   2 
and  the  T  -X  value  for  sound  speed. 

5 .   Measurement  5 

For  this  measurement,  the  projector  depth  was  8.5  cm 

and  the  hydrophone  depth  was  3.5  cm.   The  aluminum  plate 

(0.1  cm  thick)  was  buried  at  a  depth  of  5.0  ±  0.1  cm.   The 

speed  of  sound  measurements  are  summarized  in  Table  V  and 

2   2 
the  T-X  data  and  the  T  -X  data  are  shown  in  Fig.  4.11. 

Samples  of  the  received  arrivals  at  two  ranges  are  shown  in 

Figures  4.12  and  4.13. 


TABLE  V 
Speed  of  Sound  (Measurement  5) 

WAVE  SPEED  OF  SOUND  (m/s)    CORRELATION  (%)    ERROR  (%) 
SAND            1577                   99.9 

BOTTOM  1468  100  1.4 

DIRECT  1488  100  2.8 

The  parabolic  form  for  the  reflections  from  the  top 
of  the  plate  is 


45 


2 

S 
w 
K 

So 
< 


o 


X 


0                                          0 

o                             a 

<*>                               B 

*** 

WHZ 
L30< 

"""  II    II 

a  o 

- 

■    o               a 

9                                         V] 

- 

0000*2 


oooo»z 


000023 


OOO00J 

£ ( 33S0£1 


00000 1 


oooooi        oooom 


< 


o 
cu 

I 


% 


% 


o 

a 
c  o  < 

UJ  =  C0 

"J  tl  II 

a 

-  ew 

a  c 

. 

■» 


on 


ooo 


o«»  o»*  o»» 

OZSfOl 


00  » 


OK 


2       „2 


Figure  4.11.   T-X  and  T  -X   Plots  (Measurement  5) 
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Figure  4.12.   Measurement  5  at  X  =  28.2  cm, 


Figure  4.13.   Measurement  5  at  X  =  4  8.8  cm. 
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T   =   0.0450  X2  -  0.213  X  +  407 


and  the  form  for  the  reflections  from  the  top  of  the  sand 
is 


T   =   0.0463  X2  +  0.562  X  +  339 


The  incident  angle  for  determining  the  ray  parameter  was 

—  fi 

33.4  degrees  so  that  the  ray  parameter  P  was  3.70  x 10 

The  ray  parameter  analysis  gave  9.59  cm  for  AX  and 

71.5  ys  for  AT;  therefore,  the  sand  layer  depth  was  com- 
puted to  be  4.58  cm—about  six  percent  below  the  actual 
value.   The  calculated  value  for  the  sound  speed  in  the  sand 
layer  was  1904  m/s  which  was  significantly  higher  than  any 
other  of  the  determinations. 

B.   OVERALL  RESULTS 

In  all,  seven  experiments  were  performed  but  the  data 
from  two  of  the  experiments  was  not  good  enough  to  analyze. 
The  other  five  experiments  have  been  reported  above  and, 
except  for  number  5,  the  results  were  quite  consistent  with 
reasonable  errors. 

Table  VI  gives  a  summary  of  the  results  of  the  five 

experiments  that  were  analyzed  and  compares  the  results  of 

2   2 

the  T  -X   method  and  the  ray  parameter  method.   The  average 

error  in  determining  the  depth  of  the  layer  was  3.5  percent 
so  this  measurement  can  be  done  effectively;  however,  the 
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TABLE  VI 
Summary  of  Brief  Results 

MEASUREMENT 

CONTENT  12        3        4        5 

U    ACTUAL  3.5(+.05)  3.3(+.l)   2.3(+.l)   1.8(+.l)  5.0(+.l) 

Q 

yg   EXPERIMENTAL      3.56     3.73     2.42     1.74     4.58 

J     ERROR  (%)         0.3      9.7      0.8      0.0      6.5 


|  T2-X2  METHOD  1662  1647  1649  1582  1577 

^  u  RAY  PARAMETER  1791  1474  1578  1647  1904 

O  w 

TjjS.  DIFFERENCE  163  154  49.6  18.9  275.9 

w  ERROR    (%)  10.0  9.5  3.0  1.2  16.9 


error  in  the  sound  speed  calculation  (for  the  sand  layer) 
was  fairly  large  (about  8  percent) .   Consequently,  the 
sound  speed  determination  is  not  quite  as  reliable. 

In  Table  VI,  the  difference  means  the  values  between  the 

2   2 
average  of  the  five  T  -X   results  and  the  ray  parameter 

method's  value. 
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V.   DISCUSSIONS  AND  CONCLUSIONS 

The  objective  of  this  research  was  to  localize  buried 
objects  (e.g.,  thin  aluminum  plates)  in  water-saturated  sand. 
In  order  to  accomplish  this  objective,  we  needed  to  predict 
the  sound  speed  and  the  depth  of  the  overlying  layer.   The 

results  are  summarized  in  Table  VII. 

2   2 
The  average  speed  of  sound  by  the  T  -X  method  is  1472 

m/sec  for  the  direct  wave,  1446  m/sec  for  the  bottom  reflected 

wave,  1623  m/sec  for  the  wave  through  the  sand. 

Comparing  the  above  sound  speeds  with  the  value  given  by 
Reference  3  (1481  m/sec  at  20  °C) ,  we  had  0.6%  error  in  the 
direct  wave  speed  and  2.3%  error  in  the  bottom  reflected 
wave  speed. 

Reference  3  gives  the  sound  speeds  1540  m/sec  for  coarse 
silt  and  1730  m/sec  for  quartz  sand.   The  sand  used  in  this 
experiment  is  #30  mesh  sand  (washed,  kilndried  and  graded) 
and  the  measured  sound  speed  is  1623  m/sec. 

We  used  linear  regression  for  finding  the  speed  of  sound 

2   2 
with  the  T  -X   method.   Since  most  correlation  values  of  the 

2   2 
linear  regression  are  over  95%,  the  T  -X   method  is  reasonable 

for  determining  the  speed  of  sound. 

By  using  the  ray  parameter  method,  we  determined  AX,  AT 

and  Az.   The  average  error  between  the  actual  layer's  depth 

and  the  computed  depth  is  3.5%.   The  error  range  is  between 
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TABLE  VII 
The  Detailed  Results  of  Total  Experiment 

MEASUREMENT  12  3  4  5 

DBTAHS  SSL(STR        3-5<±.05>   3.3(±.l)     2.3(±.l)     1.8(±.l)     5.0(±.l) 


SPEED 
OF 


98.3      94.3      95.2      99.7      99.9 
1661      1647      1649      1582      1577 


fj    SOUND  99.4      99.6      99.9      99.7      100 


CM 
X 

I 
CM 


(m/sec) 

& 
LINEAR 


B 


1416      1419      1471      1458      1468 
99.8      99.9      99.8      99.9      100 


§  CORRELA-         1463      1437      1510      1461      1488 

TION  (%) 


£ 

1 

ft 

CD 

£  3 

B    ERROR  (%)  10.0       9.5       3.0       1.2      16.9 


P  (10   Sec/m) 

3.32 

3.69 

3.00 

5.30 

3.70 

AX  (cm) 

5.47 

3.87 

2.53 

5.8 

9.59 

AT  (ysec) 

51.4 

54.0 

33.9 

40.4 

71.5 

Az  (cm) 

3.56 

3.73 

2.42 

1.74 

4.58 

ERROR  (I)1 

0.3 

9.7 

0.8 

0.0 

6.5 

2 

C-,  DIFFERENCE 

(m/sec) 

163 

154 

50 

19 

276 

#  Average  error  of  measuring  Az  =  3.5% 

#  Average  error  of  c,  difference  =8.1% 


This  error  is  the  value  when  comparing  the  experimental  Az  with  actual 
layer's  depth. 

2  2  2 

c,  difference  is  the  difference  between  c-,  of  T  -X  and  the  value  of 

the  ray  parameter  method. 

This  error  is  the  value  when  comparing  the  c-,  difference  with  the 
average  value  of  T2-X^. 
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0%  and  9.7%.   The  ray  parameter  method  is  one  of  the  most 

useful  for  finding  an  unknown  layer's  depth.   But,  the  average 

2   2 

error  of  c1  between  the  T  -X   and  the  ray  parameter  methods 

is  relatively  high  (8.1%). 

In  every  measurement  the  direct  wave's  speed  of  sound 
is  always  faster  than  that  of  the  bottom  reflected  wave. 
These  values  should  be  the  same  because  both  are  measures  of 
the  sound  speed  in  water.   This  discrepancy  is  probably  a 
result  of  the  difficulty  in  determining  the  true  horizontal 
distance  between  the  acoustic  center  of  the  projector  and 

the  hydrophone  for  different  angles  of  the  projector.   In 

2   2 
some  cases,  the  difference  in  c,  between  the  T  -X   and  the 

ray  parameter  methods  is  over  10%.   This  may  result  from  the 
choice  of  curve-fit  for  the  T-X  data.   We  tried  several  ways 
to  find  the  best  curve  fit,  but  the  computed  results  were 
quite  different  from  the  expected  ones  except  for  those  from 
parabolic  fit. 

We  had  some  limitations  on  the  measurements  because  of 
the  use  of  a  small  tank.   In  particular,  there  were  inter- 
fering reflections.   So,  some  of  the  signals  arrive  at  the 
same  time.   That  made  the  analysis  of  the  layer  speed  of 
sound  and  depth  very  difficult.   Also,  the  pulse  length  was 
too  long  to  be  able  to  resolve  the  thin  layer.   Projector 
ringing  limits  the  pulse  length  so  that  it  cannot  be  made 
shorter  than  a  certain  duration.   Because  of  these  problems, 
this  experiment  would  be  easier  to  do  in  a  larger  body  of 
water . 
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If  we  do  this  experiment  in  the  ocean  instead  of  small 
water  tank,  then  identifying  the  wave  signals  is  much 
easier  because  of  less  multipath  interference.   In  the  ocean, 
we  can  try  ASPER  (airgun  sonobuoy  precision  echo  recorder) 
developed  by  the  Hawaii  Institute  of  Geophysics  as  one  of 
the  recommended  methods  [Ref .  1] . 
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